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ABSTRACT 


The electric and magnetic fields associated with the lightning return stroke are expressed as 
a convolution of the current waveform shape and the fields generated by a moving charge of ampli- 
tude one (i.e., the Lienard-Wiechert solution for a unit charge). The representation can be used to 
compute the fields produced by a current waveform of non-uniform velocity that propagates along 
a filament of arbitrary, but finite, curvature. To study numerically the effects of linear charge 
acceleration and channel curvature two simple channel models are used: the linear and the hyper- 
bolic. 
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APPLICATION OF THE LIENARD-WIECHERT SOLUTION TO A LIGHTNING RETURN 

STROKE MODEL 


by Robert Meneghini 
INTRODUCTION 

The purpose of the paper is to provide an alternative description of the lightning return 
stroke and the electric and magnetic fields associated with it. The point of view adopted is that the 
charge and current densities can be expressed as a collection of moving and stationary charged parti- 
cles, where the position (and therefore the velocity and acceleration) of each is assumed to be given 
as a function of time. For the lightning return stroke, charge and current distributions of this kind 
result in fields that can be interpreted as the fields produced by a moving charge of amplitude one 
(i.e, the Lienard-Wiechert field solution), multiplied by the current waveform and integrated over 
what might be termed an effective delay time. For typical return stroke parameters, the sources of 
the radiation are shown to depend primarily on the channel curvature and the rapid acceleration of 
charges near the ends of the channel. Since the formulation can be used to analyze channels of arbi- 
trary, but finite, curvature and non-uniform current waveform velocities, it appears to be somewhat 
more general than either the moment [Hill, 1968, 1969; Mamey and Shanmugam, 1971 ] or the 
piece-wise linear [LeVine and Meneghini, 1978a, 1978b] field representation. 

The organization of the paper is as follows. First, the charge and current densities are ex- 
pressed in terms of a collection of discrete charged particles. By a limiting process, continuous ver- 
sions of these sources are obtained from which the scalar and vector potentials and the fields are de- 
rived. To complete the solution, the return stroke is modeled simply as a transfer of charge between 
two points [Uman et al., 1975] , The boundary conditions appropriate to this model then determine 
the fields. The radiation portion of the fields can be integrated in the case of a finite linear filament 
if certain restrictions are placed on the velocity of the charges. Although the results are derived in 
an approximate manner, they compare well with previous work and serve as a guide to the integral 
representation for the field as the appropriate limit is approached. 

In the flnai section of the paper, numerical results are presented for the electric field as func- 
tions of the waveform velocity and the channel curvature. For this purpose two channel geometries 
are used: the linear and the hyperbolic. 
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Fields Produced by Charge Motion Along a Filament 

The charge density, p, and current density, j, associated with a particle of charge Q 0 can be 
described by the equations [Jones, 1964] 

p(x',t) - Q 0 5(x'-f(t)) 
j(x',t) - Q 0 f'(t) 5 (x' — f(t)) 

where 5 is the Dirac delta function, x' is an arbitrary position vector and f(t), f'(t) are, respectively, 
the location and velocity of the charge as a function of the time t. For a collection of n+1 charges 
(Q 0 , . , . , Q n ) where the position of the (j+l)th particle at time T is the same as the position of the 
jth particle at time T — At, the charge and current densities become 

p(x', t)» s Qj$(x'-f(t -jAr)) (1) 

j=0 

j(x', t) * 2 Qjf'(t-jAr)6(x'~f(t-jAr)) (2) 

j“0 J 

Dividing and multiplying each term of (1) and (2) by At and letting n -*■<», At •+ 0 and Qj -* 0 in 
such a way that Qj/AT is finite, then 

p(x',t) =J l(T)5(x'-f(t-T))dT (3) 

j(x', t) t(r) f'(t — t) 5 (x' — f(t — t)) dr (4) 

where t has units of amperes and is defined by 

t(r) * lim Qj/At (5) 

At-* 0 
Qj-tO 

The electric E and magnetic H fields that result from the sources p and j can be determined by means 
of the potentials <j> and A. where 

E = _V0 - 3A/9t 
H .-L 

Mo 


1 


VxA 


( 6 ) 

(7) 
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with 


4rr« 0 *(x,t) « f f |d5l2 S(t'-(t- 
x' t# 

4jre 0 c 2 A(x, t) * 1 1 §(t' - (t - r/c)) 

x' t 


r/c)) dt' dx' 


dt'dx' 


( 8 ) 


(9) 


The integrals in (8) and (9) extend over all space and time. The quantity r is the distance from the 
source point, x', to the observer, x. Substituting (3) and (4) into (8) and (9) yields 


4rre o 0(x 
4rre 0 c 2 A(x 


M ' j. 


t (r) 5 (x* — f(t* - r)) 5 (t' — (t — r/c) ) drdt , dx , 
x' t' r r 


( 10 ) 


,».///■ 


t (r )f (t* - r ) 5 (x' - f(f - r )) 6 (t' - (t - r/c) ) ir (u) 


In these and the following equations the limits on r and t' are taken to be ( ~ <» ,®°). To complete 
the solution, the boundary conditions are then imposed on the fields themselves. 

Assuming that the orders of integration in (10) and (11) can be interchanged, the volume 
integral can be done directly, giving [Jackson, 1962] 


4rre, 


4?re o 0(x, t ) = / ~ (t ~, t M dt' dr 

r t r 

,c2A(x,t) * If j(iIl(<Lzsm^SLzSM 

' - r- 


dt' dr 


where 


with 


r = x — x 


x' * f(t' - r) 


( 12 ) 

(13) 

(14) 

(15) 


Substituting (12) and (13) into (6) and (7) and carrying out the differentiations yields (Appendix) 

E(x, t) = L- J t(r)( Fi + Fo)dr 
4 ir<- * “ 


( 16 ) 


o r 
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and i r 

H(x, t) - — / i(r) r .< ( F i + F,) dr 
4/r Jr 

(17) 

where 


Fj ■ (r -v/c) (1 - (v/c) 2 )/p'-'r 2 

(18) 

*2 • r x ( (r - v/c) x v' ) /p 3 c 2 r 

(19) 

The distance r and the corresponding unit vector r are given by 


r - |x - x'l 

(20) 

r a (x - x')/r 

(21) 

with 


x' = f(t - |x - x'|/c - r) 

(22) 

In addition, the velocity v, the acceleration v', and the quantity p are 


v = f'(t - r/c - r) 

(23) 

v' = f"(t - r/c -r) 

(24) 

p = 1 _ v . r/c 

(25) 


where the prime(s) on f denote differentiation with respect to the argument u * t - v/c - r, 


To compute the fields from (16) and (17) for a given observer x and time t, the integration 
is performed over all r. Since the functional form of f is known and if we assume that it is single 
valued and defined for all values of its argument then at each r a corresponding value of x' can be 
found from (22). This immediately yields r; the velocity and acceleration are then obtained from 
(23) and (24). 

To complete the solution for the Helds, the following boundary conditions are imposed: 


f(u) = x' b 

u < 0 

(26) 

f(u) = x' t 

u>t m 

(27) 

< 

‘S' 

w 

It 

O 

u < 0 and u > t m 

(28) 

v'(u) = 0 

u < 0 and u > t m 

(29) 
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and 

c(u) = 0 u < 0 (30) 

t(u) **Q u -> o° (31) 

where x'^ , x' t are the distances from the origin to the “bottom” and “top” of the channel and 

u * t - r/c - r (32) 


To apply these boundary conditions, we use equations (22) and (26) which imply that whenever 
t > t - r/c, then x' a xy But this, in turn, gives r * |x - xy s r^, which is the distance from 
the bottom of the channel to the observer. Therefore we conclude that if r > t - r^/c then 
x' * xy From this relationship and (28) and (29) it follows that v(t - r^/c - r) = v'(t - rb/ f(! ~r) = 0 
when r > t ~ r^/c. 

Similarly, equations (22) and (27) imply that x' = x'* when r < t ^ t m - r t /c, where 
r t * fx — x' t | is the distance from the top of the channel to the observer. This condition and (28) 
and (29) give v(t - r t /c - r) * v'(t - r t /c - r) * 0 for 0 < r < t - t - r t /c. 

At those values of r for which the velocity and acceleration are zero, the integrands of (16) 
and (17) simplify considerably, Explicitly, 


(33) 

(34) 


A , 1 


Ww 

F 1 + F 2 = V r t 2 


Therefore, the fields can be written in the form 


r > t - r h /c > 0 


0 < r < t - t - r fl /c 
m t 


E(x, t) = 


| r / l fr) ( p i + F 2^ dr + ^b Q b (t ~ r b //c)/r b 2 


H(x, t) * JL 
4 rr 
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where 


r i * max(0, t - 

•m ~ V°> 

(35) 

t-, * max(Q, t — 

v°> 

(36) 


r *» 


Q b (t - r b /c> - _ 

I t(r) dr 

(37) 


t 2 


<V* " ! m " r t /'> 

C r, 

- J 5 i(r)dr 

(38) 


* / 0 


r b - lx-x' b l ; 

r b - (x-x' b )/r b 

(39) 

r t ■ |x - x' t | ; 

r t = (x - x' t )/r t 

(40) 


The model of the return stroke implied by the boundary conditions and the sources given in (3) and 
(4) is simply the transfer of charge from the point x' b to the point x' t [Uman et al., 1975] * In par- 
ticular, for times less than r b /c the observer measures a static electric field arising from a point charge 
of magnitude Q^(0) a distance r b away, where Q b (0) is the total charge of the system. From the 
perspective of the observer, and for times greater than r b /c, charge is extracted from the channel 
end and begins to propagate along the filament. For times exceeding t + r t /c, charge accumulates 
at the channel top with a magnitude given by (38). As time tends to infinity the entire charge of 
the system becomes concentrated at the top end point of the channel and the observer measures 
a static electric field given by r t Q t (<*)/r t “, where Q t («») = Q b (0). 

Models of this kind have been used in a number of studies [Uman et al., 1975; McLain and 
Uman, 1977] . Modifications of the model have been discussed by Price and Pierce (1977) and Lin 
et al. (1980). Despite deficiencies, the model appears adequate for gaining insight into the induction 
and radiation portions of the field and their dependence on observer location, channel structure, and 
the shape of the current waveform. Equations (33) and (34) are free space field solutions. To 
account for the presence of a perfectly conducting ground plane at z = 0, the following image con- 
tributions are added to (33) and (34) respectively, to give the total field: 
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* Ht(x, t)«--£- 
4 n 

~ ~ A A 

The vectors Fj , F, can be defined by ( 1 8) and ( 1 9) if r, r, v and v' are replaced by rj, rj, Vy, and v'j 
respectively, where 


r I * 1* - X'i! J rj * (x - x' r )/rj 

(43) 

x'j * (1 -* 2 z z) • f(t - r/c ~r ) 

(44) 

A A 

Vj * (1 - 2 z z) . f'(t - r/c - r) 

(45) 

v'l 3 (1 - 2 Z Z) • f"(t - r/c - r) 

(46) 

A A A A A A 

1 » xx + yy + zz 

(47) 


/ 

r l 


f A A ~ 

- i(r) r, x (F 


+ F^)dr 


(42) 


the unit vectors rj b , rj t are directed from the “bottom” and “top” of the image filament toward the 
observer, i,e., 


where 


A 

r Ib 

A 

r It 



= (x - 

- x', b )/|x - x 'it>l 

(48) 

= (x " 

■X', t )/|x - x’ It | 

(49) 


A A 


“ x b 1 

-2 zz, x' It * x' t - 2 zz . 



In the numerical calculations to be presented, the observer is taken to be at the surface of perfectly 
conducting ground plane. The total fields (actual plus image) are therefore 

E t (x, t) - 2 z z • E(x, t) (50) 

H t (x, t) = 2(1 — z z) • H(x, t) (51) 

where the fields on the right hand side of these equations are given by (33) and (34). 
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Special Cases 


There exist several cases for which the integrals in (33) and (34) can be evaluated. To obtain 
the Lienard-Wiechert solution for a single charge that moves along the path f(t), the current is written 
as i(r) * Q 5(r). Equations (33) and (34) can be integrated to give 


4 ne 0 E(x, t) 


qW 

Q(F, + F 2 ) 

<*W 


t < r b /c 

r b /0 < ‘ < t m + r t'° 
■ > t m + r ,/c 


and 


H(x, t) 


Q rx(F[ + F,) r b /c < t < t m + r,/c 
0 Otherwise 


(52) 


(53) 


where F| and F 2 are given by (18) and (19) if in (20) through (25) r is set to zero. For example, 
the r appearing in F j , F-, and (53) is given by 


r = |x - x'| 


where 


x' * f(t - |x - x'|/c) 


Since the velocity and acceleration of the charges are assumed to be zero at the ends of the filament 
and continuous everywhere, the fields are continuous as well. Equations (52) and (53) are the 
Lienard-Wiechert field solutions for a charge, initially stationary with respect to the observer, which 
traverses a finite distance before again coming to rest. 


To investigate the behavior of the radiation fields (i.e. the portion of the E and H field 
associated with the net power out of a volume enclosing the sources), we obtain from (33) and (34) 


"rad 


4tre 0 c- 


,/ 


Ti 


H 

p3 r 


dr 


H 


rad 


1 

4trc 


r 7 -> 

/ 


i(r) rxh 

p3 r 


dr 


(54) 

(55) 
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with , , 

h ■ rx j (r - v/c) x v'| (56) 

For the purpose of identifying the sources of the radiated field it is convenient to express h in terms 
of the channel curvature and arc length. The velocity and acceleration can be written in the forms 
[Hay, 1953; Stoker, 1969) 


v(u) * s'tj (57) 

v'(u) * s"fj + K s' 2 j 0 (58) 

where s is the arc length along the filament and s' * ds/du, s" 3 d"s/du” where u 3 t - r/c - r. The 
unit vectors jj ,tn and are- the tangent vector, principal normal and binormal to the filament, 
respectively, and form a right-handed orthonormal basis; the quantity K is the curvature. 


Substituing (57) and (58) into (56) gives 


b*Ks'< rx [(r x j^) - s'j^/c] > + s"( r x (r xj\)) 


(59) 


/Vs in the case of a single charge, two sources of the radiation field can be distinguished [Panofsky 
and Phillips, 1962; Jones, 1964] . If the speed of the charges is a constant v , then s' = v Q and 

t • 

s" 58 0 so that 


h 



A 

r x 


[(rxt 2 )-(v 0 /c)f 3 ]j 


(60) 


i.e. the integrands of (54) and (55) are proportional to the channel curvature. On the other hand, if 
the channel is linear, then K = 0 and 


h = s'"r x (r xtj) (61) 

For this case, the integrands are proportional to the rate of change in the speed with respect to the 
time t - r/c — r. Finally, if channel is linear and the speed constant then the radiated fields are 
zero. 
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For the case of a linear filament the integrals in (54) and (55) can be evaluated in an approx- 

A 

imate manner if certain restrictions are placed on the speed of the charges. Letting ft be a constant 
vector along the direction of propagation, then the velocity appearing in the integrals of (54) and 
(55) can be written as 

v(u) « ft v(u) (62) 

A 

where u * t - r/c ~ r. Since v'(u) * ft v'(u), then from (56) and (62) 

h m r x (r x ft) v'(u) (63) 


Making the change of variable w * v(u) and using the fact that 

_ dr/d t * r j_v/p_ (64) 

where r is defined by (20) and (22) and p by (25), Wien 


v'(u)dr = - p dw 


(65) 


If (63) and (65) are substituted into (54) a simpler integrand results. In general, the integral cannot 
be evaluated because land r are functions of w. If it is assumed, however, that at the bottom end 
point of the filament the charges are accelerated from rest to a constant velocity v Q over a distance 
much smaller than r b and further that u » 0 (and therefore i(r) ** i(t - r b /c)) over this same distance 
then (54) can be written approximately as 


: rad 


(x, t) = 


A 

- r b * 


(r b *$) 


0 V 


t(t - r b /c) 


/* v ° * 

J (>- r t 


2 w/c)"~ dw 


( 66 ) 


At the top of the channel a similar contribution is made by the deceleration of charges from v Q to 
rest over a distance small in comparison with r { . This contribution and (66) can be integrated to give 
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Similarly, 


*W X - <> * - J” |‘<* - f b / c> 'b x ® / Pb r b - ,(t - r t/ c - t m ) 't Pt r t| 

(68) 

where 

- 

P b » 1 -(v 0 /c)r b • t 

(69) 

P t *l - (v 0 /c)r t • 2 

(70) 


and t m is the channel length divided by v Q . 

The fields given by (67) and (68) are appropriate to a current pulse, having a rectangular 
velocity versus time curve, which propagates along a linear filament. Results similar to these and in 
the latter reference identical to them, have been derived by several different approaches [Uman, 
1975; Uman et ah, 1978; Le Vine and Meneghini, 1978a] . The electric field signature predicted by 
(67) is a familiar one. For times greater than r^/c but less than t m + r t /c, the magnitude of E fa( j 
is proportional to the current waveform c(t - r^/c) with a polarization given by r b X (r b X 2). For 
times greater than t m + r t /c, each component of the field is a weighted sum of the current wave- 
form with a time delayed and inverted image of itself. It is worth noting that as v Q approaches c 
and if 1 — ^ 0, 1 - ? t • 2 0, then the induction terms of (33) and (34) tend to zero and 

previously derived expressions for the total electric and magnetic fields are recovered* [Le Vine and 
Mmeghini, 1978a]. 

Equations (67) and (68) can be used to approximate the radiation fields for a current wave- 
form traversing a piece-wise linear channel comprised of n elements of lengths Lj , . . .,L n and direc- 
tions 2 , where the speed of propagation of the current waveform is 

V(t) = »„[U(t) - Wt - 5 Lj/v 0 )] 
and where U is the unit step function. 

The electric field, for example, is simply the sum of individual terms of the form of (67) 

[Le Vine and Meneghini, 1978b] and can be written 


*The static terms of the two solutions can be matched if the same boundary conditions are used. 
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-r b / c )r b x(r b *$j)/r b p b 



n—1 i(t — 

+ V - 

rj/o— Tj) A 

i 

<-T!> 

X 

+ 

l 



j*l 

*1 * 
r j 

Lp«J,J>1) 

pG‘,j)J 


i(t - r t /c - 
r. d. 

- j) w,** n >} 


(71) 


where rj is the distance from the jth channel transition point to the observer and 


- i L Jv 


k*l 


k' o 


P(iJ) * 1 — (v 0 /c) (rj • £j) 


(72) 

(73) 


Note that r t is the distance from the top of the nth channel segment to the observer, while r b is the 
distance from the bottom of the first segment to the observer. 


Terms in the summation of (71) correspond to radiation arising from the channel transition 
points. It is not obvious that these contributions are identical to the limiting form that would be 
obtained from (54). By “limiting form” we mean the field produced under conditions of a current 
pulse of constant speed and a channel geometry in which the curvature is non-zero only near iso- 
lated transition points. Unfortunately, we have not been able to carry out the integration of (54) in 
this limit. Numerical comparisons between (54) and (71) are discussed in a subsequent section 
of this paper. 


Channel Geometries 

a) Channel tortuosity*, hyperbolic channel 

One of the simplest continuous and differentiable channel models that can be used to study 
the effects of curvature is the hyperbola. In addition to its analytic simplicity, it also lends itself to 
a comparison of a two element piece-wise linear channel model if the elements are chosen to coin- 
cide with the asymptotes of the hyperbola. In fact, if the equations for the asymptotes are given 
then there exists a family of hyperbolas having these asymptotes which includes the degenerate case 


ORIGINAL PASS £9 

OF POOR QUALITY 


of the piece-wise linear channel, i.e. the asymptotes themselves. Since the curvature of the channel 
is the change in angle with respect to arc length, the integrated curvature of any hyperbola depends 
only on the angle subtended by the asymptotes*. This property provides a convenient way to quan- 
tify the effects of curvature on the radiated field. 

Assuming that the hyperbola lies in the x-z plane with the transverse axis along the x axis 
and with a origin displaced along the vertical (z) direction, then 

x2 _ < z _ d> 2 , (74) 

a 2 b2 

In the subsequent discussion and in the numerical results, we choose the 0 < z < 2d segment of 
branch of the hyperbola for which x > 0. 


The square of the speed, v, of a charged particle moving along the hyperbola is 

v~(t) = (dx/dt)' + (dz/dt)“ (75) 


To find the trajectory of the charges (or, equivalently, the trajectory of the current waveform) as a 
function of time, (74) is differentiated with respect to t. Squaring the result and using (74) and (75) 
to eliminate x" and (dx/dt)" gives 


kz-ir+ t 2 \ i/2 wo 

(u-d)2 + b2/ (a7 + bV/2 


dt 


(76) 


where 


7 = b~/(a" + b - ) 1 ' - 


(77) 


Integrating (76) from z(0) to z(t) yields [Gradshteyn and Ryzhik, 1965! 

L 1 


•1 


r 

* / o 


v(t)dt = t 0 + Cl sgn ( z “d)< — F(f3 1>q )-b E(0 p q) 


+ |z - d| 


b : + (z-d)3 \ 1/2 

7 * + (z — d)“/ 


(78) 


♦Since the channel model uses only a finite portion of the hyperbola, this statement is not strictly 
correct. For the parameters used in the numerical calculations, however, it is a good approximation. 
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where F, E are elliptic integrals of the first and second kind, respectively, 

i 

* Ffo.k) * f* (1 ~ k 2 sin 2 xr 1/2 dx (79) 

J o 

E(0,k) ■ /** ( 1 - k“ sin 2 x) 

J o 

The constant time, t , is given by 

‘o ■ c l(-7 F« 2 ,q)-bE0J 2 ,q)4 

with 

Cj = (a 2 + b 2 ) 1//2 /v Q b 
q * a/(a 2 + b 2 )^ 2 
#1 * Tarf^fld - z| /b) 

/?2 = Tan” 1 * (d/b) 

a nd sgn(x) =| 1 x > 0 

(-1 x < 0 

The constant v Q , defined below, has dimensions of velocity. 

To complete the specification of charge motion along the channel, the speed and acceleration 
are assigned so that they are zero at the channel end points and continuous everywhere. Choosing 
the speed to be a cubic polynomial in t both for 0 < t < and t 2 < t < t m and an exponential 
fortj < t < t 2 provides a sufficient number of constants to satisfy these requirements. Note 
that t m is the time taken for a charge to transverse the enti channel. One reason for choosing this 
particular model is that it allows us to obtain approximations to the velocity versus time curves com- 
puted by Leise and Taylor ( 1 977). 

The result is 



(82) 

(83) 

(84) 

(85) 

( 86 ) 
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v(t) 


v 0 (t/tj) 2 (aj a 2 (t/tj)) 
- U(t-tj) 


0 < t < tj 

t ] <t<t 2 (87) 

v o ((t “ V /‘f> 2 (a 3 + a 4 (t " t m )/t f ) e ~ ^ * 1 } t 2 < t < t f 


v o e 


m 


where v Q = v(t j ) is the maximum speed and 


'f * ‘m- 1 : 


aj = 3 + 0 tj ; a 2 - - 2 - 0 t| 

®3 ~ 3 — = 2 — 


( 88 ) 

(89) 

(90) 


Expressions for the z components of the velocity and acceleration can now be found from (76) and 
(87) 



dz/dt 

d'z/dt" 


bv(t) / (z -d ) 2 + b 2 \ 1/2 

(a 2 + b2) l/2 ^( Z _ d) 2+ 7 2 j 

v'(t)v z (bv(t )) 2 (z - d) ( 7 2 - b 2 ) 

" ■ 1 — - + - " ■ ■■ 

v(t) (a- + b^) ((z - d )2 + y -)2 


(91) 


(92) 


where v'(t) is determined by differentiating (87). Equations for x, v x , v r are derived from (74) and 
(75) by a similar procedure. Since x is more conveniently computed from (76). the results are given 
here only for v x , v' x which are 


v 


X 


av(t) 

(.2 + b 2)‘/2 



m 

sgn (z - d) 


(93) 


v'(t)v x (a v(t )) 2 (a 2 - a 2 )x 


y' = - 4 . 

v X T 


V(t) 


(a- + b”) (x“ — 0 “ )“ 


with 


a - aV + b 2 )*/ 2 


(94) 

(95) 
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The foregoing equations can now be used in (33) and (34) to generate the fields for arbitrary obser- 
vation points and times. Since the evaluation of the integrals is somewhat complicated and since 
simplifications are possible, it is worth indicating how the calculations are made. 

To find z as a function of t, equally spaced values of z are taken from 0 to 2d and the corre- 
sponding value of t is obtained from (78). To express z in terms of equally spaced time intervals, a 
linear interpolation is performed. The x(t) value is then found from (74) with x > 0. Once x(t) and 
z(t) are known, the velocity and acceleration are computed from equations (87) and (91) through 
(94). These quantities then determine f (t), f'(t), f"(t). 

As indicated earlier, for each sample point r in the numerical integration of (33) and (34), r 
can be found from the equations 


r =» jx-x'l (20) 

with 

x' « f(t - |x -x'|/c - r ) (22) 

An easier way to evaluate r is possible, however. Letting 

x' = f(T) (96) 

then (22) can be written 

f(T) = f(t — |x — f (T)|/c — r ) ( 97) 

The function f is assumed to be one to one so that the arguments of f can be equated: 

T - t - |x-f (T)|/c - r (98) 

Since t,x and fare known and r is fixed, (98) can be solved numerically for T. From T,x' is found 
from (96) which determines, in turn, r, v and v' , 

b) Linear Channel Model 

For the case of a linear channel, the specification of f(t), V (t) and f"(t) is straightforward. 
Letting (Xj, 0, Zj), (x-», 0, z 0 ) be the bottom and top end points of the channel then the equation 
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for any point (x, 0, z) on the channel is 

z - Zj * s(x - Xj) (99) 

where 

s * ( z 2 z 1>/( X 2 “ x l > » x 2 * x l ( 100 ) 

Solving for x and z as a function of the time t yields 

x(t) *Xj+(l+s 2 ) f 0 v(t)dt ( l01 ) 

z(t) a z, + s(l +s 2 ) ” 1/2 f v(t)dt (102) 

These equations determine f(t) * (x(t), 0, z(t)); f (t) andf"(t) are then obtained by successive differ- 
entiations. Note that if x^ 35 ^ z ( { ) can be obtained from (102) by letting s -►<». As in the hyper- 
bolic channel, the speed v(t) is assumed to be given by (87). 

Numerical Examples 

The cases presented illustrate some of the characteristics of the field solutions given by (33) 
and (34). In all examples, the observer is located on the surface of a perfectly conducting plane; 
thus, only the normal electric field and tangential magnetic field are non-zero, The current waveform 
is assumed to be 


i(r) - 30,000 (e“‘ 04r - e“ 2r ) (103) 

where r is expressed in microseconds and t (r) in amperes. In all figures the quantities v , tj and t^ 
in (87), (88) are taken such that v Q * c/3 and t j = t^ Each curve consists of 400 equally spaced 
samples in the time domain. The CPU time required for a representative curve is approximately 40 s 
on an IBM 3081. 

To show the effects of changes in speed independent of the channel curvature, a linear chan- 
nel model is used. Examples are shown in figures 1 and 2. For both figures the channel is vertical 
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with bottom and top end points atx*y*z*0andx*y*0, z*4 km, respectively. The observer 
is located at x * 100 km, and y * z * 0. In figure I the curves A, B, C correspond to values of t j , as 
defined in (87), equal to 16, 8, 2 ms respectively. In all cases 0 * 0, Curve D was generated by means 
of (67); all others from (33). As mentioned earlier, curve D represents the field in the limiting case 
t j 0 and is proportional to a sum of two current waveforms where one is displaced in tinvs and in- 
verted with respect to the other. Increasing t j has two effects: the curves are ‘stretched' and the 
peak amplitude of the field is reduced such that the time integral of the magnitude of the field re- 
mains approximately constant. It should be noted that in the case of a linear filament (67) has been 
shown to approximate (33) well as long as tj is less than about. 1 ms. This remains true for different 
channel orientations and propagation speeds, v Q . 

In figure 1 the parameter 0 as defined in (87) equals zero, In figure 2, tj is fixed at 4 MS 
while 0 is varied. The curves A, B and C correspond to values of 0 equal to 2 x 10^, 1 .5 x 10^ and 
10^ s”* , respectively. The speed versus time curves for these cases are shown in the inset. Since 
tj and the maximum velocity, v Q , are constant, the peak field occurs at approximately 4 ms in all 
cases. 


Increasing 0 tends to narrow the width of the waveform and to decrease the positive over- 
shoot. This latter effect is caused by a decrease in acceleration near the channel top; for example, 
at 0 * 10^ s~ l the charges must be decelerated from v = 6.3 x 10^ ms - ^ to zero within 4 ms while 
fo:r 0 * 2 x 10^ s -1 they are decelerated from the smaller value v = 2.5 x lO' 7 ms“^ to zero within 
the same time span. 

To explain the narrowing of the initial portion of the field, note that for 0 7 = 0, charge accel- 
eration at the channel bottom produces a negative z directed field. For 0 =?= 0, the charges, having 
reached their maximum velocity, are then decelerated, i.e., v'(t) = - which results in a posi- 

tively directed E,,. The total field, therefore, is the sum of negative and positive z components, the 
latter contribution beginning 4 ms after the first. The result is that the width of the negative portion 
of the waveform narrows as 0 increases. It should be noted that the field signatures given by curves 
A and B are qualitatively similar to the radiated waveform proposed by Leise and Taylor (1977). 
This is not surprising since we have essentially proceeded in a reverse direction by starting with a 
velocity curve qualitatively similar to their derived curve and then using this to compute the field. 

To study the effects of tortuosity on the fields, the speed of the current waveform is chosen 
to be a constant, v Q , in figures 3 and 4: i.e., t j = 0, t^= 0 and 0 * 0 in (87), This situation is not the 


$ame as obtained by taking the limit tj * tf*+ 0 which yields radiation fields of the form given by 
(67), Rather, the assumption is simply used to eliminate the sources of radiation from the ends of 
the filament and isolate curvature effects. 

In figure 3 are shown three curves of the z component of the electric field for an observer at 
x * 100 km, y * z * 0. The channel geometry is a hyperbola, given by (76), with x>0, 0<z<3km 
(d * 1.5 km) and a * b. Curves A, B and C correspond to values of ‘a* equal to 500, 100, and 20 m, 
respectively. The channel geometries are shown in the inset. Since the ratio b/a is constant the 
asymptotes of the channels are identical and therefore the total integrated curvature of each of the 
channels is approximately constant. As ‘a’ decreases the field is seen to tend toward zero for times 
less than 25 ms After 25 ms the shape of the curve C corresponding to a * 20 m is similar to the cur- 
rent waveform. This behavior can be qualitatively explained by noting that as ‘a’ decreases, the 
only portion of the channel where significant charge acceleration occurs is near the ‘bend’, i.e,, the 
apex of the hyperbola. In the limit of small ‘a’ the radiation field, (which is the dominant contrib- 
utor to the total field at this distance), originates entirely from this point and is proportional to the 
current waveform, t. Several examples were run for ‘a* less than 20 m, down to a minimum of 2 m. 
The corresponding field change, however, is small. 

Figure 4 shows the effect of changing the slope of the asymptotes. The channel geometry is 
again given by (74) with x>0, 0<z<3km but now with a = 2b. Curves A, B. and C correspond 
to values of ‘a’ equal to 500, 100, and 20 m. Notice that amplitudes of the field are increased as 
compared to those of figure 3. This is reasonable in the sense that the integrated curvature of the 
channel has increased over that assumed in figure 3. Since the integrand of (54) is proportional 
to curvature for the v « v Q case then, in general, larger fields can be expected. 

The combined effects of channel curvature and linear acceleration on the electric and mag- 
netic fields are shown in figures 5 and 6. For both figures the channel geometry is hyperbolic with 
x > 0, 0 < z < 3 km and a = b = 100 m. The velocity is again given by (87) withj3 3 0 and t j = 2 ms. 
Curves A, B and C correspond to observer distances of 1 km, 10 km and 100 km along the x-axis. 
For curves A and B in figure 5, the static contribution is evident. This is a consequence of the 
simple charge transport model used and represents the accumulation of charge at the top end of 
the channel. Since the bottom of the channel is on perfectly conducting plane, the static con- 
tribution from this endpoint is zero. For the magnetic field, plotted in figure 6, the static contri- 
bution is, as shown by (34), zero. 
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A number of numerical comparisons were made between the fields predicted by the piece- 
wise linear representation of (71) and the Lienard-Wiechert solution of (33), The comparisons 
were made m the far field and therefore (33) and (54) are essentially the same, To make the 
geometries as close as possible n, in (71 ), was set equal to 2 and for (33) a hyperbolic channel 
was chosen with parameters V and ‘b’ as small as possible without incurring sampling errors. In 
addition tj * t f in (87) was chosen to be small to approximate the delta function acceleration at 
the end points inherent in (67) and (71), As indicated earlier, the radiation arising from the end 
points of the channel are in good agreement for the cases considered. On the other hand, the fields 
arising from the channel transition points showed differences which generally increase with increas- 
ing channel curvature. For example, when the curvature is small, 2a = b, the two solutions are in 
good agreement. At a =* 2b, however, the ratio of the Lienard-Wiechert to the piece-wise linear 
solution is approximately two for v 0 between c/10 and c/ 1 *25. As noted before, we have been un- 
able to analytically evaluate (54) in the limiting case of a piece-wise linear channel. Some qualita- 
tive understanding of the differences between (54) and (71) can be gained, however, by noting that 
in deriving (71) the velocity and acceleration were assumed to be colinear and therefore the term 
hj « r x (v/e x v') in the integrand of (54) was zero. On the other hand, as we let the parameters 
‘a’ and V in the equation for the hyperbola approach zero, the geometry will tend toward a piece- 
wise linear channel. In general, however, h j will be non-zero and therefore will contribute to the 
Lienard-Wiechert solution. 


Discussion and Conclusions 

A form of the Lienard-Wiechert solution has been used to determine the electric and mag- 
netic fields associated with a commonly used lightning return stroke model, An apparent advantage 
of this form of solution is that parameters of the current waveform can be assigned independently 
of the channel geometry, Because of this, a class of geometries can be studied in the context of the 
same return stroke model. For example, one type of geometry that can be analyzed is a cubic spline 
approximation to photographed lightning channels. 

Effects of linear acceleration and curvature were studied by means of two simple channel 
geometries: the linear and the hyperbolic. For the special case of a linear filament and a rectangu- 
lar velocity versus time profile, the solution is in good agreement with previous work. Significant 
changes in the signature appear, however, when the shape of the velocity profile is altered. 
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In the case of the hyperbolic filament, the radiated Fields are functions of the curvature, ob- 
server location and current waveform shape, Only in cases where the non-zero curvature is concen- 
trated about a small segment of the channel do the signatures approximate the current waveform 
shape, Preliminary numerical results also have shown that if the velocity profile of the current wave- 
form is rectangular and the non-zero curvature occur only near isolated points along the channel, 
then the Lienard-Wiechert and the piece-wise linear representations are in fair agreement in many 
cases. In particular, the two solutions give nearly identical results for the radiation arising from the 
channel end points. As to the radiation from the channel transition points, discrepancies between 
the two solutions are observed that become more pronounced with increasing channel curvature. 
One of the reasons for the differences appears to be the assumption used in deriving the piece-wise 
linear solution; that the velocity and acceleration are colinear. 

It should be mentioned in conclusion that several other issues pertaining to the Lienard- 
Wiechert representation have not been answered: whether the Field given by (33) is equivalent to 
the Panofsky and Phillips solution as applied to the return stroke model by Leise and Taylor (1977); 
and whether a form of the Lienard-Wiechert representation can be derived for more complicated 
models of the return stroke (Price and Pierce, 1977 ; Lin et al., 19801 . 
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APPENDIX - CALCULATION OF THE FIELDS 


The electric E ?,nd magnetic H fields are given by 


E * - V * — 3A/3t 

H - JL V*A 


where <p and A are, from (15) and (16), 


with 


4ffe 0 _ _ 

0 T t' 


// t(r)5(t' - (t - r/c)) dt# 


dr 


4rre 0 c-A 


If i(r) f'(t' - r) 6 (f - (t - r/e» dt' dr 


r t' 


(A!) 

(A2) 


(A3) 


(A4) 


r* I x — f (t' —r) 


The limits of integration on t' and r are (— <*, »). 


(A5) 


Taking the gradient of (A3) and interchanging the order of integration and differentiation yields 

47re o V0 * J c(r) ( l^(r) + lg(r)) dr (A6) 

where 


I A (r) 


--/ 4 

t' r 


5 (t' - ft -r/c)) dt' 


l B (r) = j -j- Vd (t' - (t - r/c)) dt' 
t' * 

By an application of the chain rule T75 ft' - (t - r/c)) can be written as 


(A 7) 
(A8) 
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r 5 '(t' - (t - r/c))/c where the prime on the delta function denotes differentiation with respect to 
its argument. 

Following Jackson (1962), the change of variable w * t' + r/c is made in (A7) and (AS). 
From (A5) we find that 


dw/di' * 1 - v • r/c * p 


(A9) 


with 


v * f'(t' - r) 


Equations (A7) and (A8) become 


I 


* “/ 


r 5 (w - t) 


dw 


pr 


... f U' 

® f) ' J ~ 

l^(t) can be integrated immediately to give 


(w-t) 


prc 


dw 


l A (r) ■ -(r/pr 2 ) c 


(AiO) 


(All) 


(A 12) 


(A13) 


where the subscript o means that all variables within the brackets are to be evaluated at t' 3 t- r/c. 
Explicitly, 


with 


r = |x - x'l 


p = 1 - v • r/c 


x' = f(t - |x -x'| / c - r) 


v = f'(t - |x — x'J/c - r) 


(A 14) 
(A15) 
(A16) 
(A17) 
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Integrating (A 12) by parts gives 


'B 


(r) 


!- (— il 

(dw Vprc/jo 


Using the facts that 


JL (IV p - 1 - (- 

dw \pr / dt' \pr y 


-^(pr)” 1 * -Cpr)“ 21 

dt' 


( v 2 a rir • v'\ 

v . r 

c c ) 


* A 

dr 

d? 


r~ l (v - r(v . r)) 


then 




U 2 


cl 


r (v"/c - v * r - r r • v'/ c) + (v •» r(v • r))p 


A,.. A, 


J JO 


From (A6), (A14) and (A22), 

i , r 

r (v-/c ~v • r - r r • v'/c) + p(v - r(v r)) 
The calculation of 3 A/3 1 follows closely that for V<p. The result is 


4ire 0 Vp =» f i(r) { — i- 

U 3 r 2 c 


4jre 0 c“ 3A/3t =ji (r) | prv' ~ v v 2 /c - v •r-rr^v'/c | 


dr 


Inserting Vtf into 3A/dt into (Al) and using the identities 


A 

r x 


(r - v/c) x v' 


= r • v' (r - v/c) - v' p 


( r - v/c) (1 - (v/c) 2 ) - r p 2 - p(v - r (v • r))/c - (v 2 /c - v • r) (r - v/c)/c 


then 


(A 1 8) 

(A19) 

(A20) 

(A21) 

(A22) 

dr (A23) 

( A24) 

(A25) 

(A26) 
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1 

E(x,t) * ~~ Jt.( r)(F,+F,)df 

(A27) 

where 

F 1 * I -7- (r-v/c) (1 -(v/c) 2 | 

I p 3 r 2 ) 0 

(A28) 


F, - <-rr r x ((r-v/c)x v'] 1 
1 p 3 c 2 r ) 0 

(A29) 

To obtain the magnetic field, the curl of (A4) is taken which gives 



4tre 0 c 2 Vx A(x, t) =* -fi (r) (l c (r) - l d (r)) dr 

(A30) 

where 

» c (r) =/ f ^ ~ r)Xf 5'(t'-(t-r/c))dt' 
J rc 

(A31) 


f f'(t' - r) X £ 

= J r 2 5(t'-(t -r/c))dt' 

(A32) 

Proceeding as before and using (A2) and the identities 



-r x | r x [(f-v/c)x v'] j = rx j (r • v')v/c + pv' | 

(A33) 


-c 2 £x [ (“r — v/c) (1 — v“/c“)] = (vxr)(v*-c*) 

(A34) 

then 

H(x, t) = — - ft (r) <|r x (F 1 + F 7 ) j dr 
4 ; t J 1 - )o 

(A35) 


To simplify the notation in the text the subscripts o have been omitted. 
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Figure 4. Same as figure 3 except that a = 2b. 
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